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^ ■ Abstract 
00 '• 

\^ , In this letter the exphcit form of evolution operator of the Tavis-Cummings 

' model with three and four atoms is given. This is an important progress in quantum 

■ optics or mathematical physics. 

O ' The purpose of this paper is to give an explicit form to the evolution operator of 

^ ' Tavis-Cummings model (^J) with some atoms. This model is a very important one in 
Quantum Optics and has been studied widely, see [2] as general textbooks in quantum 
^ '. optics. 

\ We are studying a quantum computation and therefore want to study the model from 



this point of view, namely a quantum computation based on atoms of laser-cooled and 
trapped linearly in a cavity. We must in this model construct the controlled NOT gate 
or other controlled unitary gates to perform a quantum computation, see |^ as a general 
introduction to this subject. 
5^ I For that aim we need the explicit form of evolution operator of the models with one, 

two, three and four atoms (at least). As to the model of one atom or two atoms it is 
more or less known (see jlj), while as to the case of three and four atoms it has not been 
given as far as we know. Since we succeeded in finding the explicit form for three and 
four atoms cases we report it (0)- 
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The Tavis-Cummings model (with n-atoms) that we treat in this paper can be written 
as follows (we set ^ = 1 for simplicity). 

An n 

^ i=i 1=1 



where ui is the frequency of radiation field, A the energy difference of two level atoms, 
a and are annihilation and creation operators of the field, and g a coupling constant, 

'i 1 



and L = 2". Here cr,-^'*, cr,- ^ and of^ are given as 



cjf ^ = I2 ® • • ■ ® I2 ® ® I2 ® • • ■ ® I2 (« - position) G M(L, C) (2) 
where s is +, — and 3 respectively and 

^+=(0 0)' ^-=(1 0)' ^=^=(0 -1)' ^^=(0 l)- 
Here let us rewrite the hamiltonian (^. If we set 

s^=T.^t\ 5-=E^^^ s,=^Y.^?\ (4) 

1=1 i=l ^ i=l 

then (|TJ can be written as 

E = ujIl ® a)a + AS'3 ® 1 + g (S+ ® a + S_ ® a)^ =Hq + V, (5) 

which is very clear. We note that {S'+, S^, S^} satisfy the s'u(2)-relation 

[53, 5+] = 5+, [53, SJ\ = -5_, [5+, 5_] = 2^3. (6) 

However, the representation p defined by p(cr+) = 5+, p(o"_) = S^, p{a^l2) = 5*3 is a 
reducible representation of 5^(2). 

We would like to solve the Schrodinger equation 

ijU = HU= {Ho + V) U, (7) 

where ?7 is a unitary operator (called the evolution operator). We can solve this equation 
by using the method of constant variation. The result is well-known to be 

under the resonance condition A = u, where we have dropped the constant unitary 
operator for simplicity. Therefore we have only to calculate the term (jH)) explicitly, which 
is however a very hard task ^. In the following we set 

A„ = 5+ O a + S*. O a"^ (9) 



^The situation is very similar to that of the paper quant-ph/0312060| in [J] 

2 



for simplicity. We can determine e '^^^^ for n = 1 (one atom case), n = 2 (two atoms 
case), n = 3 (three atoms case) and n = 4 (four atoms case) completely. 



One Atom Case In this case A in is written as 



By making use of the simple relation 
with the number operator we have 



a 
at 



Bl/2 



(10) 





^ 




f N + 1 \ 


[ 






V N 



-itgAi 



^ (-1)" 

( 



2n+l A2n+1 



COS [tg^jN + 1 



V 



.sm tgVAf , 



o(2^ + l)' 
.sin(tgv^]v+T) 

cos {tg^/N^ 



{tgr^' A? 



(12) 



We obtained the explicit form of solution. However, this form is more or less well-known, 
see for example the second book in [2^ . 



Two Atoms Case In this case A in 



is written as 



A, 



f a a \ 

at a 

at a 

V at at y 



(13) 



Our method is to reduce the 4 x 4-matrix A2 in ()13|1 to a 3 x 3-matrix Bi in the 
following to make our calculation easier. For that aim we prepare the following matrix 



T 



then it is easy to see 

Tt^aT 



/ 



V 



/ 1 \ 
V 001; 



V2a 
V2at V2a 
V2at J 



where -Bi = J+ ® a + J_ ® at and {J+, J_} are just generators of (spin one) irreducible 
representation of Q. We note that this means a well-known decomposition of spin 
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Therefore to calculate e ^^^"^^ we have only to do e . Noting the relation 



(2{N + l) 2a2 
2(2A^+1) 
V 2(at)2 2N 

2{2N + ?,) 

2{2N + 1) 

2(2iV- 1) 



fii = DB,, 



and so 



= D"-'Bf for n > 1, 



2n+l 



for n>0 



we obtain by making use of the Taylor expansion 



-itgBi 



1 + E 



-1)" 



-1)" 



{tgr+'Bl^+' 



nt^l (2^)! ^ " ' n=o(2n+l)! 

/ 1 + i±i/(iV + 1) -^h{N + l)a 2ife/(iV + l)a' 
-ih{N)a^ l + 2f{N) -ih{N)a 



where 



fiN) 



-1 + cos (tg^2{2N + l) 



h{N) 



sm 



{tg^2i2N+l) 
V2N+1 



A, 



CcLSG 
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in © 


is 


written as 
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(14) 



(15) 



(16) 



We would like to look for the explicit form of solution like (fT2|) or p5|l . If we set 
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then it is not difficult to see 
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at 
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V3at 





2a 












2at 





v^a 
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v^at 





/ 



1/2 



1/2 



-B3/2 / 



This means a decomposition of spin \® \ ® \ = |©|©|. Therefore we have only to 
calculate e"**^^^/^^ which is however not easy. In this case there is no simple relation like 
(fTT|) or ((HI), so we must find another one. 

Let us state the key lemma for that. Noting 



-^3/2 



r3 

-D3/2 



2y3a2 



\ 

2 



2V3a 



/ 3A^ + 3 

7N + 4: 

2v^(at)2 7A^ + 3 

\ 2v^(at)2 3N J 

( V3(7A^ + ll)a 



V 



v^(7A^ + 4)at 


6(at)3 



and the relations 



R2n+1 
^3/2 





20iVat 




-83/2-6372; 



20(iV+ l)a 


\/3(7A^ - 4)at 



rj2n+2 _ r2 R2n 
-^3/2 ~ -D3/2-D3/25 



6a^ 


v^(7iV + 3)a 




we can obtain B'^^ s-^d B^^^ like 



B 



,2n 
3/2 







f «„(iV + 2) 

/3„(iV + l) 

2v^en(iV)(a^)' 

V 2VSUN-l){a^y 



2y/3^r^{N + 2)a 


7n(iV) 





2 



5n(iV-l) 



R2n+1 
-°3/2 



where 



/ V^Pn{N + 2)a 

v^/3„(iV+l)at 2e„+i(iV + l)a 

2e„+i(iV)at 

V 6UN - l)(at)3 V3^n{N - l)at 



QUN + 2)a^ \ 


y37n(A^)a 
J 



J 

(17) 

(18) 



an{N) = {v+Xl - t;_A!!)/(2v^), /3„(iV) = (w+A*^ - i^_A^)/(2v^), 
7„(iV) = (t;+A''! - t;^A'[)/(2v^), <5„(iV) = (w;+A!! - w;_A^)/(2v^), 
UN) = iXl-X-_)/i2Vd), 



and A± = X±{N), v± = v±{N), w± =w±{N), d = d{N) defined by 

X±{N) =5N± ^Jd{N), v±{N) = -2N - 3 ± ^Jd{N), w±{N) = 2N - 3 ± yJd{N), 
d{N) = IQN'^ + 9. 

Then by making use of ()17p and ()18|) we liave 



oo / 1 \ n. 



n=0 (2^;^ 

/ /2(iV + 2) 



3/2 



2n+l j^2n+l 



oo / 1 \n 

-V3iFi(A^ + 2)a 2y3/ii(A^ + 2)a2 -6iifo(A^ + 2)0^ \ 
/i(iV + l) -2ii/i(iV + l)a 2y3/ii(A^+ l)a2 

2^3/11 (A^) (at) 2 -2iHi{N)a) fo{N) -^/3iFo{N)a 

\ -6tHo{N - l){a^Y 2V3hi{N - l){a^y -V3tFo{N - l)a^ f_i{N-l) j 

(19) 



-V3iFi{N + l)i 



wfiere 



/2(iV) 

fo{N) 
f-iiN) 

Fi{N) 

Fo(iV) 
Hi{N) 
Ho{N) 



v+iN)cositg^KiN)) - v^{N)cos{tg^X^{N))j /i2^d{N)), 
w4N)cos{tg^KiN)) - ^_(iV)cos(t(7^A_(iV))} /{2^d(N)), 
v+{N)cos{tg^\^{N)) - v^{N)cos{tg^K{N))\ /{2^d(N)), 



w+{N)cos{tg^X_{N)) - w^{N)cos{tg^^ \+{N)) \ /{2^d(N)) 



cos{tgsJX+{N)) - cos{tg^X4N)) \ /{2^d{N)) 



^^Id^sm{tg./X4N)) - ^^l£Lsin(t^7yrM) | /(SV^)- 



sin(t^?JA_(iV)) 



V-iN) . 



, , - , , sin(t<7v'A+(iV)) /(2yrf(iV)), 

A_(iV) ^A4iV) J 

^X+{N)sin{tg^X+{N)) - ^ X.{N)sm{tg^ X^{N))\ /{2^d(N)), 



^ sm{tg^X+{N)) - sm{tg^X^{N)) \ /{2^d{N)). 

^X+{N) y^A_(iV) ' 



6 



Four Atoms Case In this case A4 in (jH)) is written as 
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If we set T as 
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then it is not difficult to see 
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This means a well-known decomposition of spin |©|©|®| = 0©1©0©1©1©2. 

Since we have calculated e~^^3Bi (|P) we have only to do e~**^^^, which is of 

course hard. The proof is similar to that of the three atoms case, so we state only the 
result. 

exp {-itgB2) = 



( f2{N + 2) h{N + 2)a'^ fco(A^ + 2)a^ \ 

/i(A^ + l) ho{N + l)a^ 

hi{N){a^y foiN) h^i{N)a^ 

ho{N-l){a^y /_i(iV-l) 

V A;o(iV-2)(at)4 /i_i(iV - 2)(at)2 f^2iN-2)J 
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(21) 



where 



UN) = 1 + 4(iV - l){K/A+)(cost(7yA+ - 1) - {u^/X^){costg^X^ - 1)}/Vd 
fi{N) = (m+ cost(7-y/A+ — n_ costg\J~xl) / \fd, 

fo{N) = l + 2{{v+w+/\+){co^tg^jK-l) - (t;_t/7_/A_)(cost(7yAr - l)}/v^, 
f-i{N) = {u+ cos tg\J~\l — U- cos tg\J~X^) / Vd, 



/_2(iV) = 1 + 4(iV + 2){(n+/A_)(cost(7v/Ar- 1) - {u_ / \+){costg^+ - I)} / V~d, 
h{N) = 2{{v^/X+){costg^-l) - iv_/\_){costg^ - 1)}/Vd, 
hQ^N) = (cos t g — cost g^JJ^) 

h_,{N) = 2{(«;+/A+)(cost(7^- 1) - {w^/\^){costg^ - 1)}/Vd, 
fco(iV) = A{{l/X+){costg^-l) - {l/X_){costg^-l)}/Vd, 
and 

Fi(iV) = {{u+/^^)smtg^-{u^/^) smtg^}/Vd, 
F_i(iV) = {{u+f^) smtg./L-{u-/^) smtg^}/vQ, 
Hi{N) = 2{{v+l^+) sintg^-iv^/^) smtg^}/Vd, 
Ho{N) = {(l/^)sint(7V^- (I/^aT) smtg^}/Vd, 
H^,{N) = 2{{w+/^) smtg^-{w^/^) smtg^}/Vd 
, and d = d{N), \± = \±{N), u± = u±{N), v± = v±{N) and w± = w±{N) defined by 

\±{N) = lOAT + 5 ± 3^fd(N), u±{N) = ^(-3 ± ^/d(N)), 

v±{N) = ^(2iV - 1 ± ^/W)), ^±(iV) = ^(2iV + 3 ± v/d(iV)), 
d{N) = AN^ + 4A^ + 9. 



A comment is in order. We note that in the process of calculation we used Math- 
ematica to the fullest (a calculation by hand might be "painful"). We would like to 
generalize the results in this paper to the cases of more than four atoms. However, it is 
not easy to perform a calculation due to some severe technical reasons. There is a (big ?) 
gap between the four atoms and the five ones. 



We obtained the explicit form of evolution operator of the Tavis-Cummings model for 
three and four atoms cases. This is a big progress in quantum optics or mathematical 
physics ^. Therefore, many applications to quantum physics or mathematical physics will 
be expected, see for example papers in [4 . 

We can also apply the result (s) to a quantum computation based on atoms of laser- 
cooled and trapped linearly in a cavity [Hj. 

We conclude this paper by making a comment. The Tavis-Cummings model is based 
on (only) two energy levels of atoms. However, an atom has in general infinitely many 
energy levels, so it is natural to use this possibility. We are also studying a quantum 
computation based on multi-level systems of atoms (a qudit theory) jjj. Therefore we 
would like to extend the Tavis-Cummings model based on two-levels to a model based 
on multi-levels. This is a very challenging task. 

^To obtain an explicit solution for some (intersting) model is still important 
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